人年の計算

library(survival)

## Loading required package: splines

data(leukemia)
# time → 追跡期間に対応する変数（標準では日数を想定） scale :
# 人年の単位（標準では日を年に変換するためscale=365.25） offtable :
# どのカテゴリにも属さなかった人年の和（エラーチェック用）
res <- pyears(time ~ x, leukemia, scale = 1)
res

## $call
## pyears(formula = time ~ x, data = leukemia, scale = 1)
## 
## $pyears
##    Maintained Nonmaintained 
##           423           255 
## 
## $n
##    Maintained Nonmaintained 
##            11            12 
## 
## $offtable
## [1] 0
## 
## attr(,"class")
## [1] "pyears"

# Surv()を使えばイベント発生数もカウント data.frame=T
# とすれば，データフレームで出力
res <- pyears(Surv(time, status) ~ x, leukemia, scale = 1, data.frame = T)
res

## $call
## pyears(formula = Surv(time, status) ~ x, data = leukemia, scale = 1, 
##     data.frame = T)
## 
## $data
##               x pyears  n event
## 1    Maintained    423 11     7
## 2 Nonmaintained    255 12    11
## 
## $offtable
## [1] 0
## 
## attr(,"class")
## [1] "pyears"